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We measure the time dependence of the scintillator light-emission pulses in Nal(Tl) crystals at 
different temperatures, after activation by gamma rays. We confirm that there are two main non- 
exponential components to the time decay and find that their amplitude ratio shows Arrhenius 
temperature dependence. We explain these nonexponential components as arising from two com- 
peting mechanisms of carrier transport to the Tl activation levels. The total light output of the 
Nal(Tl) detectors shows a linear temperature dependence explained by our model. 



PACS numbers: 

I. INTRODUCTION 

Doped alkali-halide scintillation crystals, developed 50 
years ago, continue to be the workhorses of outdoor 
gamma spectroscopy because of their size, cost, and per- 
formance. Recently, there has been a growing interest in 
Nal(Tl) detectors for use in Homeland Security applica- 
tions [TJ [2] - However, a well-known, problematic feature 
of Nal(Tl) detectors is the temperature-dependence of 
their light yield. The core of this problem stems from 
the temperature dependence of the shape and amplitude 
of the light pulses emitted from the scintillator for a given 
energy of the incident ionizing particle. Although multi- 
ple exponential components describing the shape of the 
light pulse in time have been reported [3], the common 
understanding for a Nal(Tl) light pulse assumes a single 
dominant exponential component with a temperature- 
dependent decay-time constant [H [5j |6j |7j [8], based on 
the effective three-level activator model described in [5]. 
This is inconsistent with the well known fact that be- 
low 60°C, the current pulses cannot be fitted well with 
a single exponent, even at a fixed temperature. Until 
now, it has been also generally accepted that the tem- 
perature behavior of the total light output of Nal(Tl) 
is nonlinear and has a broad maximum below room 
temperature. (3JIH [HE]. 

Recent experimental results obtained by using a novel 
method and device demonstrated that this experimen- 
tally observed nonlinear temperature behavior arises 
because of a convolution between the temperature- 
dependent shape of the light pulse and the pulse re- 
sponse of the shaping circuitry [TU]. It was shown ex- 
perimentally that when a gated integrator is used, the 
whole light output shows only linear temperature depen- 
dence over a wide temperature range (-30°C to +60°C) 



p~0l flTl H2l IT3] . In the same work, two main compo- 
nents of the shape of the light pulses, with a tempera- 
ture redistribution between their amplitudes, were found. 
Importantly, it was shown in jTTJ Q2] that the slow 
component is negligible above room temperature, but it 
produces up to 40 percent of the total light at — 20°C 
and lasts several microseconds. It was also shown that 
at high temperatures, only one almost exponential decay 
component exists, consistent with [6]. 

We measured the time dependence of Nal(Tl) cur- 
rent pulses at different temperatures. We show that 
rather than one or more exponential decay components, 
there are two dominant nonexponential light compo- 
nents, corresponding to two competing mechanisms of 
carrier transport to the Tl activation levels. We model 
the transport and reproduce the experimental data with 
nonlinear rate equations. The data may also be fitted 
with two exponential decay components using fixed fast 
and slow time constants and an Arrhenius temperature- 
dependent redistribution between the amplitudes, but 
the accuracy is approximately 5 percent. 



II. EXPERIMENTAL RESULTS 

We placed a standard 2-in x 2-in. Bicron Nal(Tl) de- 
tector in an environmental chamber. The chamber tem- 
perature changed at a rate of 2°C/hr, and each set of 
measurements at a given temperature was taken after an 
8-hr hold time to allow good thermal equilibration. We 
digitized the photomultiplier tube (PMT) current pulses 
with a 14-bit digital scope, manufactured by Gage Ap- 
plied Sciences Inc., connected via a 50 f2 resistor directly 
to the anode of the PMT for each temperature. The 
temperature coefficient of the PMT was almost constant 
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over the investigated temperature interval [2] and did 
not change the shape of the light pulse. The area of 
the pulses was selected to correspond to an excitation 
around 662 keV. We normalized all pulses to the unit 
area to study the effect on pulse shape alone. In Fig. 1, 
we show the time dependence of the emitted light out- 
put on a log scale, where the two components are re- 
solved and the temperature-dependent redistribution be- 
tween the amplitudes of the slow and fast components 
is readily seen. Note that the emission spectrum of this 
doped inorganic scintillator shows that the emitted light 
is at a single wavelength, characteristic of a transition 
between a single excited state of the activator and its 
ground state. Thus the two components of the light pulse 
do not represent decays from different energy levels, but 
rather, they represent two different pathways for the sec- 
ondary electrons and holes to the activator levels. We 
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FIG. 1: Current pulses measured at 

-20°C,-10 C,0 C, 10 C , ,20 C,30 C,50 C , . Each one 
normalized to its maximum value so the second component 
may be clearly resolved at lower temperatures. 

normalized the area of the current pulses to remove the 
temperature-dependence arising from a competition be- 
tween the scintillation transitions and the nonradiativc 
phonon-assisted transitions. This competition take place 
at the excited Tl level, where the electrons and holes 
recombine [3]. The nonradiative, phonon-assisted tran- 
sitions decrease the total amount of the collected light, 
but they play this role after the secondary carriers, elec- 
trons and holes, are recombined at the (Tl + )* level. The 
nonradiative, phonon-assisted transitions do not signifi- 
cantly change the time shape of the current pulse. This 
can readily be understood at very higher temperatures 
where only one exponential component of the current 
pulses exists [6] and the time shape of the pulse is almost 
temperature independent. On the other hand the non- 
radiative, phonon-assisted transitions are responsible for 



the strong temperature dependence of the peak position, 
and this dependence can be observed at arbitrarily higher 
temperatures. At such high temperatures, the shape of 
the current pulse does not change with temperature, but 
the area of the pulse and the corresponding peak position 
for a given energy of the incident particle becomes less 
and less with increasing temperature.. 

In our measurements of the light pulses from Nal(Tl), 
we found an Arrhenius dependence of the ratio between 
the amplitudes of the two dominant components with 
which we fit the data. We found a similar dependence 
in CsI(Tl) using the numerical data published in [IS], as 
shown in Fig. [2] We found that 

— oce ( (1) 

where Qi(T) and Q2(T) are the amplitudes of the fast 
and the slow components, Qi(T) + Q2{T) — constant, 
and AE is a phenomenological activation energy between 
the STE level and the excited Tl level. 
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FIG. 2: Activation dependence of the total amplitude of the 
current pulses for Nal(Tl) (squares) and Csl (Tl) (triangles), 
on a logarithmic scale. Numerical values for the amplitudes 
Qi and Q? of GsI(Tl) by courtesy of [T5]. 



III. THEORETICAL MODEL 

Because we were sure that the observed components of 
the light pulses did not correspond to impurities or dif- 
ferent energy levels, we sought to determine the primary 
transport mechanisms for secondary electrons and holes. 

In alkali-halide crystals, the electron- hole interaction 
is very strong [l6j . This leads to the formation of ex- 
citons and, because of the highly polarizable nature of 
the crystals, to a fast creation of self-trapped excitons 
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(STEs) at the very beginning of the process. This is why 
the earliest theoretical models of the transport in Nal(Tl) 
included only the motion of the excitons [T7] . There is 
considerable evidence of long-lived STEs in alkali-halidc 
crystals at low temperatures [TS]. In Nal(Tl), a redis- 
tribution between the emitted light from the STE levels 
and the emitted light from the Tl activation centers above 
T= 140 K has been observed. This means that after the 
decay of an STE, the created electron and hole usually 
reach the activator centers [19]. The suggested trans- 
port mechanism is by phonon-assisted exciton hopping 
from the STE level to the Tl level. In Nal(Tl), this type 
of transport was observed with a temperature-dependent 
rate of 10 12 - lO 8 ^" 1 ] [19]. 

On the other hand, a previously performed experiment 
on CsI(Tl) showed that the predicted dependence of the 
shape of the light pulse on the activator concentration, 
for different energies of an incident particle, is not ob- 
servable [T7| HH]- This disagreement with the existing 
theory led to the incorporation of a second kind of motion 
in the crystal: namely, the binary diffusion of electrons 
and holes [2T]. The binary diffusion in these materials 
leads to the fast creation of localized electrons and holes 
in the vicinity of the activator centers and, after that, to 
complex mechanisms of recombination. Optical measure- 
ments [25J |2S] proved that the lifetime of holes trapped 
at Tl ++ is surprisingly short: approximately 0.350 fisec 
in Nal(Tl). Therefore, we considered another possible 
mechanism of STE decay, leading to the creation of sep- 
arated electrons and holes and to the diffusion-transport 
mechanism. The mechanism we considered is the previ- 
ously discovered multiphonon-assisted dissociation of the 
STE [24]. The corresponding rate is again temperature 
dependent, and its value is in the same range as the rate 
of the hopping transport. 

The simple relationship shown in Eq.([T|) suggests that 
there are two transport mechanisms in competition for 
reaching the activated site. We propose that one is a 
phonon-assisted hopping transport of the STE to the Tl 
level, and the other is a multiphonon-assisted dissocia- 
tion of the STE followed by single-carrier transport of the 
electron and hole to the Tl level, as shown in Fig[3] In the 
first channel, the STE reaches the activator levels (Tl + ) 
via phonon-assisted hopping, making (Tl + )* excited lev- 
els that creates the fast component of the light output. 
The decay time of this component would be mainly the 
lifetime of the (Tl + )* excited levels. We can infer that 
the combination of the lifetime of the STE level and the 
STE transport to the activation center is relatively fast. 
Note that the decay rate is determined also by the non- 
radiative, phonon-assisted transition between (Tl + )* and 
the ground state (Tl + ). 

We associate the second (slower) component with bi- 
nary transport 52]. The mechanism is by STEs that ther- 
mally dissociate into electrons and holes. The electrons 
are rapidly captured at (Tl + ) levels, making (T7°) levels, 
and the holes are quickly captured at (Tl + ) levels, mak- 
ing (Tl ++ ) levels. The holes may reach the (Tl°) levels 



via diffusion and recombine with the electrons, thus cre- 
ating (Tl + )* excited levels that decay optically, as before. 
The lifetime of this process is slower than the hopping 
transport because it is a combination of the relatively 
long lifetime of the trapped hole, the diffusion time, and 
the lifetime of the (Tl + )* excited level. The two processes 
are schematically depicted in Fig. [3] 
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FIG. 3: Diagram representing the two main competing pro- 
cesses in alkali-halide crystals. Intermediate processes creat- 
ing Tl ++ and Tl° are omitted. The quenching levels in Tl + 
are not shown because the nonradiative transitions change 
only the integrated light output. They are responsible for the 
linear temperature dependence of the total light output, but 
they do not change the temporal shape of the light emission. 

Generally speaking, it is clear that both parts of the 
response are not simply exponential, being a convolution 
of the lifetimes of the trapped hole, diffusive transport, 
and the (Tl + )* level in the slow component, and of the 
lifetime of the STE level, hopping transport, and lifetime 
of the (Tl + )* level in the fast component. 

We now compare this physical picture with the tradi- 
tional one that is currently used HI H]. Making an 
approximation, the established model is based on a sin- 
gle temperature-dependent decay-time constant [B] . This 
model actually describes the lifetime of the activator level 
|3J. It is supported from a simple consideration of the 
activator's center consisting of three energy levels: the 
ground level G, the scintillation level S, and the quench- 
ing level Q |9J . The calculations of the quantum efficiency 
q in this three levels-model lead to 

fc s exp(-^s-) 

q "fc s exp(-^) + fc Q exp(-M^)' 



Three levels model takes into account that at temper- 
ature T, the relative populations of the levels arc as fol- 
lows: level S ~ exp(-f£); level Q ~ exp (- Ws ^ Q ). 
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Here kg is the probability of radiative transition from ra- 
diative level S to the ground level G, and ItQ is the proba- 
bility for nonradiative transition from the quenching level 
Q to the ground level (so called internal quenching). Ws 
is the energy difference between scintillation level S and 
ground level G, and Wq is the energy difference between 
quenching level Q and scintillation level S. The quantum 
efficiency calculated in this way is a monotonically in- 
creasing function with decreasing temperature, and one 
would expect a monotonic temperature dependence of 
the total light output. 

However, experimentally, the temperature dependence 
of the light output from the Nal(Tl) scintillator was 
found to be strongly nonlinear [3j H] , possessing a broad 
maximum below room temperature, in disagreement with 
the monotonic behavior of the quantum efficiency of the 
model. An additional function f(T) (so called exter- 
nal quenching) that represents the probability of the oc- 
cupation of the activator's centers or, in other words, 
describes the temperature dependence of total charge 
Q total — Qtotai(T), was introduced to multiply the quan- 
tum efficiency in order to explain this nonlinear behavior 
[3J |3] . This function has been introduced to describe the 
nonlinear shape of the temperature dependence of the 
total light output. But as we pointed out, the nonlinear 
temperature dependence does not exist if one collects the 
total light using a gated integrator [ID]. The nonlinear 
behavior of the peak position with temperature is differ- 
ent for different shapers [TUJ , because the long component 
produces a different fraction of the whole light at differ- 
ent temperatures [TU1 ED EH E2] • Using different shapers, 
we do not gather the whole light, but we cut off part of 
the signal produced by the long component, and there- 
fore, we cut off part of the light that was produced via 
the slower pathway to the Tl level. This cutting is tem- 
perature dependent because the amount of the light that 
we cut belongs to the long component, and the amount 
of the light belonging to the long component is tempera- 
ture dependent [TUl ED E2J El] • The existing models for 
the temperature behavior of the Nal(Tl) do not take into 
account the complex shape of the individual light pulse, 
connected with the redistribution between the slow and 
fast component, or the simple linear temperature depen- 
dence of total light output that was recently measured 

ma- 

In our model, the total charge created from the in- 
cident ionizing particle has two main components. As 
we discussed above, we found experimentally that their 
amplitude ratio shows an Arrhenius temperature depen- 
dence, whereas their sum is a constant. This means that 



Q hopping (T*) 
Qdiss ociation 

(T) 



QtotaiAexp ^ 
1 + Aexp|| 

Qtotal 

l + Aexpff ' 



activator levels, and A is the ratio between the proba- 
bility for STE decay via multiphonon dissociation and 
the probability for STE decay via hopping at very high 
temperature. Both components of the total charge are 
strongly temperature dependent, but their sum does not 
depend on temperature. That is why in our model we 
have f(T) = 1. As a consequence, we have temperature 
dependence of the total light output, as in the pure three- 
level activator model. The quantum efficiency calculated 
in this model, pointed out in Eq.([2]), is a monotonically 
increasing function with decreasing temperature. This 
means that in a narrow temperature interval (kT in the 
range 0.03 to 0.02 eV), it can be approximated with a lin- 
ear temperature dependence. Therefore in such a narrow 
temperature interval (T in the range -30°C to +60°C ), 
the total collected charge Qtotal depends almost linearly 
on temperature, and because of that, we will observe an 
almost linear temperature dependence of the total light 
output L(T): 

L(T) ~ -CT+ 1, 

where C is a constant, and we measure in percent L(T) 
as it was shown in riDJ ED EU E3]- At the same time, 
we have a strong temperature dependence and redistri- 
bution between the two main components of the light 
pulse, which explains the variety of the nonlinear tem- 
perature dependences of the light output when different 
shapers are used [10] . 

To check the model, we simulated the processes consid- 
ered above using rate equations that describe the popula- 
tions of the STEs, separated electrons and holes, popula- 
tions of the (Tl ++ ), (Tl°) levels, and the final population 
of the (Tl + )* level, which gave us the time dependence 
of the light output: 
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Here, Q h opp mg {T) and Qdissociation{T) correspond to 
charges that use two different pathways for reaching the 



The equations include nonlinear bimolecular terms de- 
scribing the diffusion of the holes from (Tl ++ ) to (Tl°) 
and thermoactivated transport of the electrons from 
(Tl°) to (Tl ++ ) that are responsible for afterglow. We 
numerically solved this stiff system of ordinary differen- 
tial equations using a standard MatLab program, odel5s. 
Our variables are as follows: 
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(1) STE - density of the STEs. 

(2) e~ - density of the electrons. 

(3) h + - density of the holes. 

(4) Tl° e _ - density of trapped electrons at Tl° levels. 

(5) Tl~^+ - density of trapped holes at Tl ++ levels. 

(6) Tl+* pul - population of (TI+)* excited levels. 
The parameters are also shown below: 

(1) T hop - the lifetime at the STE level, combined with 
the time an exciton needs, via hopping, to reach the 
Tl excited level. 

(2) Tdi S s - the time an STE needs to decay to electron 
and hole via multiphonon dissociation. 

(3) T t ° rap - the lifetime of an electron before it is trapped 
at Tl°. 

(4) r t ++ p - the lifetime of a hole before it is trapped at 
TI++. 

(5) T a ftgi - the time an electron trapped at the Tl° level 
needs to reach and recombine with a hole trapped at 
Tl ++ via thermoactivation. 

(6) Tdiff - the time a hole needs, after detrapping from 
the Tl ++ level, to reach and recombine with an elec- 
tron trapped at Tl° via diffusion. 



(7) r s , 



the lifetime at excited Tl level. 



(8) - the lifetime of a trapped hole at the Tl ++ 
level. 

(9) T° trap - the lifetime of a trapped electron at the Tl° 
level. 

The most important parameters we used are well 
known: the lifetime of the trapped hole is approximately 
0.350 /isec [23], and the lifetime of the (Tl+)* excited lev- 
els is approximately 0.134 fj,sec [B]. Less well known are 
the other parameters such as the multiphonon dissocia- 
tion time, phonon-assisted hopping time, etc., which we 
varied. But provided these parameters are kept within 
reasonable ranges, the overall results are insensitive to 
their values. The values we used for T = — 20°C and for 
T 



C are shown in Table I. 



TABLE I: Parameters Used in Simulations (nsec) 
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The obtained solutions fit our experimental results 
very well, as can be seen in Fig 4. This figure clearly 



shows the overall nonexponential form of both the ex- 
perimental results and the model. 
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FIG. 4: Two of the experimentally measured current pulses 
normalized by area for Nal(Tl) at -20°C (upper curve) and 
+50°C (lower curve) are shown in grey dots. The numerical 
solutions of the nonlinear rate-equation model are shown in 
black dots. The inset shows the same pulses on a 1 /xsec scale. 



IV. CONCLUSIONS 



We have presented a model of two competing trans- 
port processes that explains the existence of two tempo- 
ral components in the scintillation of Nal(Tl). While the 
linear temperature dependence of the total light output 
is due to the presence of a nonradiative decay channel 
from the activator level, a much stronger temperature 
dependence in Nal(Tl) arises because of the temperature- 
activated type dependence of the ratio between the fast 
and slow components. These two nonexponential compo- 
nents correspond to two distinct mechanisms of transport 
to the activator levels. One of the predictions made by 
the model is that the differential light output as a func- 
tion of the energy of the ionizing particle will depend on 
the temperature. Furthermore, the temperature depen- 
dence of the shape of the light pulses, for a given initial 
energy, will depend on the level of doping because it arises 
form the competition between unimolccular and bimolcc- 
ular processes. In a subsequent publication, we will con- 
sider these topics, and we will argue that the difference in 
the transition rates between Nal(Tl) and CsI(Tl) is de- 
termined by the differences between the available phonon 
density of states at the activation energy for these two 
materials. 
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